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Abstract -In this work, we show that a gate electric field, applied in the base of the field-effect 
devices, leads to inducing spin-orbit interactions (Rashba and linear Dresselhauss) and confines 
the transport electrons in a two-dimensional electron gas. On the basis of these phenomena we 
solve analytically the Pauli equation when the Rashba strength and the linear Dresselhaus one 
are equal, for a tuning value of the gate electric field £* . Using the transfer matrix approach, we 
provide a joint description of the transport by varying the bias electric field, £b- We can flip the 
spin of the incident electrons, or block the spin-down completely. The robustness of this behavior 
is proved when £* changes by £g ± S£ g . 
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Introduction. — The main consequences of the 
nanometer-scale devices, with defined boundary condi- 
tions, are the confinement of the electric carriers in one- 
dimensional potential wells and the formation of discrete 
energy spectra. For example, when an homogeneous mag- 
netic field is applied in a nanostructure, the energy levels 
are defined by the well-known Landau levels and the as- 
sociated spectrum can be manipulated by this magnetic 
field [TJ . For a gate electric field, there are two approaches 
for studying its influence in a two-dimensional electron 
gas (2DEG) spectrum: with perturbation methods 2 4 
or with Airy functions [SHZ] ■ We found the quantization of 
the transverse energy of the electron, the thickness of the 
2DEG and, hence, the induction of the spin-orbit interac- 
tions, by taking into account the exactly solvable problem 
whose solutions are linear combinations of two indepen- 
dent Bessel-(±i) functions and simulate the gate junc- 
tion as an infinity triangular potential well for an idealized 
Spin-FET, see fig. Q]a) and b). 

The manipulation of spin by electric fields in semi- 
conducting environments has generated a lot of theoreti- 
cal and experimental research aimed at developing useful 
spintronic devices and novel physical concepts [8j|9]. The 
spin transistor elucidated by Datta and Das [10] is made 
to drive a modulated spin-polarized current. For this, the 
spin precession is controlled via the Rashba spin-orbit cou- 
pling associated with the interfacial electric fields present 



in the quantum well that contains the 2DEG. The strength 
of this interaction can be tuned by the application of an 
external gate voltage [TTHM] . The linear Dresselhaus spin- 
orbit interaction (SOI) could be induced by a gate electric 
field 7 1 15 - IS] . This work shows how the gate electric field 
can tune the Rashba and linear Dresselhaus strengths and 
the interplay between them. For gate electric field value 
where the interplay is equal to 1 and its neighborhood, 
we study the spin transport of the 2DEG in the base by 
solving analytically the Pauli equation and build the cor- 
responding transfer matrix. 



The 2DEG and the SOI terms induced by a gate 
electric field. — In this section, we shall study how the 
gate electric field, £ g , works on the base. It quantizes, 
with the appropriate boundary conditions, the transverse 
energy of the electrons, it induces the spin-orbit interac- 
tions and the creation the 2DEG. In fig. [T]a), it is shown 
an idealized spin-FET, with a semiconductor base's di- 
mensions w x , w y and w z and two ferromagnetic materials 
as the source and the drain. We ignore the magnetic field 
caused by the ferromagnetic contacts. The source-to-drain 
electric field in the z direction, £{,, is weaker; hence, we ne- 
glect its influence to induce any SOI in our approach. 

In order to study the influence of the gate electric field, 
we take into account the one-dimensional Schrodinger 
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Fig. 1: a) An idealized Spin-FET with base's dimensions w x , 
w y and w z . We consider here that T is oriented along the crys- 
tallographic axis [100] . A 2DEG at the plane yz will be formed 
into the base when a gate electric field is applied and it has a 
thickness given by the average position (x). b) The gate elec- 
tric field provokes a confinement due to a triangular potential 
well V (x) = ~e£ g x for —k%1% < x < 0. A bound wave- 
function itself must vanish outside of the triangular well and 
at the endpoints x = —k^Ib = — {^n) 2 ^ Ie s and x = 0. 
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where m is the effective mass, e is the electron charge and 
E x is the transverse energy. It is well-known that the so- 
lution of this equation is given by the linear combination 

of the Bessel-(±i) functions X (x) = x?J-\ f§X^) ^ + 
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Several electrons can be trapped in the base, 
when the gate electric field is applied. In order to under- 
stand the electric field's confinement on the base, we take 
into account a triangular potential well V (x) = ~e£ g x 
formed by V (x) = oo for V (x) < —e£ g x and for x > 0, 
as shown in Fig. Q] b), this model is equivalent to the 
used ones in refs. [5HZ]- The bound states inside of this 
triangular potential can be determined if we consider the 
following boundary conditions X (x = — K x l% ^ =0 and 
X (x = 0) = 0. The most general solution of the one di- 



mensional Schrodinger equation, satisfying the first con- 
dition, is of the form X (x) = x^J\ (§X'l B- O n the 
other hand, the second boundary condition implies that 



B = 0. Hence 



X = K x l Eg and x (0) = K x l Eg Ji 
nontrivial solutions of this problem exist only if k x has a 
value such that 



root of the Bessel--| function, 



x,n = (| r n) 3 here r n is the rath- 
i. e. Ji (r n ) = 0. In 

3 v ' 

this way, the n-bound state acquires the energy E x . n = 

2 2 2l 

(l r ™) 3 2ml 1 = (f r ™) 3 e ^g l E g - Geometrically, a bound 
wave-function itself must vanish outside the triangular 

2 

well and at the endpoints x = — (| r„) 3 l Eg and x = 0, see 

2 

Fig. [T]b). Taking into account a shift — (§ r n) 3 lE g of the 
zero of the x-coordinate, we can write the boundary condi- 
tions, X (0) = X f(|fri) 2 ^ 3 Ie b ) = 0, the one-dimensional 



Schrodinger equation ([1} is reduced and has the follow- 

,2 v 

ing structure: + jt~X = 0. The average squared 

Eg 

momentum can be calculated with the previous equation 



(Pi) = 



By using the boundary 

,2/3 



conditions x = and x — (ffn)^' ^E g anci considering 
a Sturm-Liouville problem, we can calculate the average 
position (x) 



r(f-*) j 



0* 



(2) 



dx 



With the change of variable x = (§ r n£) 2 ^ 3 Ie 9 , the previ- 
ous integrals can be written as 



Jo^ /3 [ji(^o] 2 d 



(3) 



The integral of the denominator has a well-known value, 

r i 2 

\ J± {fn) , while the integral of the numerator is cal- 
culated with the help of the following indefinite integral 



when 9 = 1/3 and (3 = r n . In this way, (x) = | (| \ 
and, therefore, the average squared momentum is written 

by (pi) = | (|r„)^ (h/l Eg f - A 2DEG at the plane yz 
will be formed into the base when a gate electric field is 
applied and it has a thickness given by the average posi- 
tion (x). This thickness is numerically equivalent to the 
one given in [7] for the ground state, r%. This thickness 
decreases monotonically when £ g grows, but its magni- 
tude increases for higher values of r„; this implies that 
the 2DEG is created only by small root values r n and 
high gate electric fields. 

We consider here that both SOI are induced by 
the gate electric field £ g . If the geometry of the base 
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Gate electric field [xlO 4 eV/cm] 

Fig. 2: Interplay t as a function of the gate electric field. 
In plot a) GaAs and the roots n = 2.9026, r-2 = 6.0327 and 
rs = 9.1705. The Dresselhaus strength increases with large 
values of r n and the Rashba one does not. When j*- > 1 
the Rashba SOI dominates the behavior of the 2DEG, while 
jpr < 1 the behavior of the 2DEG is given by the Dresselhaus 

term. Plot b) shows the interplay ^ as a function of £ g for 
three conventional III-V semiconductors InSb, InAs and GaSb 
at the ground state n. The dominant SOI in the InSb is the 
Rashba term, while the GaSb has a strong Dresselhaus one. 
The SOI terms of the InAs behave similarly as the GaAs does. 



is such that w x > w z > w y » (x) for an applied 
gate electric field, it is possible that the inequalities 

(pI) » (p 2 y ) gc ^ (pI) sc > (pI) sc can bc truc ' where 

the sub-index gc means geometrical confinement and, 
for example, (Px) gc — ( 7r ^- / ' w x) 2 for the grown state. In 
the effective-mass approximation, the electron Hamilto- 
nian of a zinc-blendes structure has a spin-dependent 
p 3 coupling called the Dresselhaus term [15] , Hd = 

f£ [<?xPx (P 2 y ~ Pf) + <?yPy (pi ~ P%) +°zPz if x ~ Py)} , 

where (3 3 d is the Dresselhaus three dimensional term, pi 
and a-i with i = x, y and z the components of the momen- 
tum and the Pauli matrices respectively. For the triangu- 
lar potential well and the ^-direction, which is consider 
parallel to crystallographic direction [100], the Dressel- 
haus term is reduced to Hd = ~^¥- (?yPy — CzPz)? 2 - We 
propose that the most important term of the linear Dres- 
selhaus interaction is written by Hd = — j- {fJyPy ~ a zPz\ 

with f3* ~ f (|r„) ^ P3d/1 2 b ■ For a spin-FET where 
w x ~ (x) the linear Dresselhaus strength could be given 

by p* = fan [| (|r„) f l/l 2 Eg + OM) 2 ] ; in this work 
we neglect the geometrical confinement. consider 
only the geometrical confinement because w x << (x). 
On the other hand, the Rashba interaction [231 has the 



Table 1: Physical properties of the zinc-blende structure semi- 
conductors, according with ref. [7] 



Semiconductor 


m 
mn 


a [A*] 


P 3D [eVA 3 ] 


GaAs 


0.0670 


4.4 


26 


InSb 


0.0136 


500 


228 


InAs 


0.0239 


110 


130 


GaSb 


0.0412 


33 


187 



VzPy) 



with a* 



aeE a 



The ratio ^ denotes the 



following form Hr = \- [pyP* 
and a the material's constant. 

interplay between the Rashba strength and the linear 
Dresselhaus one. Such interplay is given by the following 
expression 



5 _ 

3 rf: 



(4) 



it is a function of the electric field and depends paramet- 
rically on the root r n and the semiconductor nature (m, 
a and /33£>, table [1] has those values). By varying the gate 
electric field, one can modulate this interplay. Figure [2^,) 
shows the interplay ^ for GaAs, graphed in the interval 

1 x 10 3 ^ < £ a < 5 x 10 4 ^ at the roots n = 2.9026, 

cm y cm 1 7 

r- 2 = 6.0327 and r 3 = 9.1705. grows monotonically 
with the gate electric field for a fixed root r n . The hori- 
zontal dotted line point up when ^ — 1. If we take into 

account the curve r-± we can find that ^ = 1 for a value 
of the gate electric field, £*, inside the considered interval, 
but for T2 and r 3 these fields are too large. In other words, 
the Dresselhaus strength increases with large values of r n 
and the Rashba docs not. When > 1 the Rashba SOI 

dominates the behavior of the 2DEG, while < 1 the 
behavior of the 2DEG is given by the Dresselhaus term. 
Figure [2j)) shows the interplay ^ asa function of £ g for 
three conventional zinc-blende semiconductors InSb, InAs 
and GaSb at the ground state r%. The curve of the in- 
terplay for InSb is over 1, while the corresponding curve 
of the GaSb is below 1. In this way, the dominant SOI 
in the InSb is the Rashba term, while the GaSb has a 
strong Dresselhaus one. The SOI terms of the InAs be- 
have similarly as the GaAs docs. Those analitical calculus 
generalizes the numerical results of ref. [7j . 

The transport of the 2DEG under induced spin- 
orbit interaction. — Into the base, several electrons 
can be trapped by the gate electric field E = £gi and they 
form a 2DEG, which is moving in the plane yz under the 
influence of the induced SOI terms. This movement is 
described by the Pauli equation with SOI 



H$ (x, y, z) = E F <& (x, y, z) 
where H = -^p 2 -e£ g x-e£bZ+V (y) + £ f- (<r y Pz - <x 2 



(5) 



(<TyPy 
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is the Hamiltonian, p 



-iftV is the 



p-3 



J. L. Cardoso and H. Hcrnandez-Saldana 



momentum, m is the effective mass, is the bias electric 
field (£f, « £ g ), Ef is the Fermi energy and a* and (i* 
are the Rashba and the Dresselhaus strengths respectively 
and V(y) is the transverse confining hard wall potential 

(y{y) = f° r < y < it an< ^ mnm t° outside this po- 
tential region). According to [24) . the time reversal opera- 
tor T is anti-linear, i. e. T(Ti — > —<Ji, and it has the follow- 
ing structure T = ia y C, with C the complex-conjugation 
operator. The Hamiltonian with SOI terms is invariant 



under time reversal, 



0. Thus, this Hamiltonian 



can not produce spontaneous spin polarization |12j . To 
get rid of the variable x, we will consider the complete 
set of stationary states of a particle in an x-dimensional 
triangular potential well |x„ (x) = Xr/ 2 Ji Xr/ 2 ^ -B„| 

with Xn = 7^ — I - (| r n) 2 ^ 3 - I n other words, we can use 
these functions to express the wave function <f?(x,y,z) in 
the form 



-y' 



4>i (y, ■ 



(6) 



where the expansion coefficients <pj (y, z) are spinors. If 
we introduce this function <&(x, y, z) in the Pauli equation 
((SJ, we have 



*+2i (or, + ^-(§r n )^ 



9jP 



2* (^ia^l; 



</»„ (x, y) = 0, 



(7) 



here ^ rp — r» 

Imetb 



a = ma*/h 2 and /3 = m/3*/h 2 . The 
degrees of freedom of the x axis are completely separated 
from those in the yz plane. This problem can be con- 
sidered then as a two-dimensional one and its solutions 
depend parametrically on r n . We choose the confine- 
ment along x to be much stronger, for high gate electric 
fields, such that only the lowest subband (given by n) 
is occupied in this direction under all operating condi- 
tions. We only look for solutions around n, other modes 
can be solved in a similar way. For clarity, the index 
n = 1 shall be suppressed in the following expressions. 
In order to get rid of the variable y, we will consider 
the complete set of stationary states of a particle in a 
y-dimcnsional infinite potential well with boundary con- 
ditions i n (—w y /2) = L n (w y /2) = 0, which are solutions of 
the differential equation 



tly- 



2i (Pa y 



acr z ) jfe +q y , r L = 0. 



It is easy to verify that t — J2 r i r = (y) <V, 

where l, L r and c r are spinors, Y r (y) = cos (k y _ r y) for 
r odd and Y r (y) = sin (k y ^ r y) for r even, k v>r = ^p, 



k 2 

™y,r 



1 + (3 2 and R y = e i ^' y y +0 " y ^ v is a rotation 
operator. If we use the complete set of functions {t r } to 
expand <f>(y,z), we get <f>(y,z) = RyE^i^fe)^^)- 
By introducing this function in the Pauli equation ([7)l. 



multiplying by R y Y s (y) and integrating on the variable y, 
we have 



T 



2i 



+ 4^ (2a 2 + AJ S 



I (2/3 2 



AJ, 



(8) 



~ 2 V J2r [Wxlr,s - (a<T y ~ /3(T Z ) J r>a ] -jf = 



here n 2 
A = c 

by I r ,s - 
2 r w v/ 2 

Wti J- 



2mE F 



k 2 - 



2.6664 



V 



P 2 , 



2 — /3 2 , and the mixing terms are given 
: ^ S-wj/2 sin2 ^ Y r{y)Ys{y)dy and J r . s = 
^ 2 cos2r]y Y s (y)Y r (y)dy. Notice that if (3 = ±a 
the previous equations are uncoupled; in ref. )20j a similar 
case is studied. However, this symmetry would be broken 
with a magnetic field pointing in any direction. 

According to fig. [2] a) for the GaAs [fig. [2] b) for the 
InAs], the value of the gate electric field £* = 2.6126 x 
10 4 g£ {£* = 2.6735 x 10 4 g£) gives the interplay | = 1. 
If we take values in the neighborhood of £ * ± 5£ g , where 
= ±0.06, it is possible to show that (3 = a (1 ± 0.02) 

a 

and the eq. (J8} can be reduced in the first approximation 



2a 2 (1± 0.02) 



-2ia [<r y + (1 ± 0.02) cr 2 ] Z 
±i0.04a£ r [VWr,- - [<r y - tr,) J,.,,] 



(9) 



dz 



0. 



The coupling terms are not significant and, therefore, they 
can be negligible and the previous coupled eqs. become 
the following system of uncoupled eqs. 



2ia (a y +a z )± 



+ k 2 + 2a' 



0. 



(10) 

For each wavenumber k z>s , there are two propagat- 
ing physical channels: one with spin-up and another 
with spin-down. The general solutions of these equa- 
tions are given by Z s (z) = A s (z) C\ + B s (z) C2 be- 
ing A s (^) = e- ialE ^ + °*K°(h-i (0), B s (z) = 



ch* (id), c. 



+ 



4a 2 ) l\ 



and C\ and Ci arc spinors. We now match the wave func- 
tion Z s and its derivative Z' s at the borders of the base 
region. In this way, we can connect the incident waves (to 
the left-hand side of the source) with the outgoing ones (to 
the right-hand of the drain). The evolution of the spin- 1/2 
state vectors is governed by the two-channel 4x4 transfer 
matrix 



M B (z) = A" 1 



A s (z) B s (z) 
A' s (z) B' s (z) 



TA 



(11) 



where 



1 1 

* K z,sl * K z.sl 

1 

ial Eb {(Ty + o z ) 1 
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Fig. 3: Transmission coefficients as functions of the bias elec- 
tric field. In the panel a) we have the transmission coefficients 
T^fi and T^-f, while in the panel b) we have the transmission 
coefficients Tj,^ and T-1-4. By varying Sb, the transmission can 
be modulated more efficiently: we can flip the spin of the in- 
cident electrons, or block the spin-down completely, and thus 
establish a spin transistor. 
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Bias electric field [eV/cm] 

Fig. 4: Transmission coefficients as functions of the bias electric 
field. In panel a) T-f.ti T4.4 and T±f axe plotted for S£ g /£g = 
+0.06, in panel b) If,-)-, and T^f are plot-ed for 8£ g /£g = 
—0.06. The general trends discussed in fig. [3] are preserved 
and, especially, one can modulate them by varying the bias 
electric field. 



A' s and B' s are just the derivatives of A s and B s with re- 
spect to z and 1 is the unit matrix. This transfer matrix 
has the symplectic structure [24] 



M, 



a 

kb*k T 



b 

ka*k T 



(12) 



with k = i<7y. Given we are ready to calculate the 
whole Spin-FET transmission amplitude 



(a f )" 



*t,t *t4 



I.I 



(13) 



where Uj (here i and j label the spin-up (t) or spin-down 
(J.) projections) is the transmission amplitude from chan- 
nel j on the source to channel i on the drain, and a is 
the (1,1) block of the spin-FET transfer matrix Mg. The 
off-diagonal terms t-j- j, and are the transmission ampli- 
tudes for processes where an odd number of spin flips have 
taken place inside the base. The corresponding transmis- 
sion coefficients are defined by T ? J = . 

We consider electrons of energy Ep — 78mcV in- 
jected into a InAs base structure, with £* = 2.6735 |^ 
for j — 1. In this report, we will keep fixed the fol- 
lowing geometrical variables: w x = 5 x 10 _1 ^im and 
w y — w z = 5 x 10~ 2 fim. We are interested here on the 
visualization of spin-transitions based on the transmission 
coefficients. We shall start considering the effect of vary- 
ing the bias electric field on the transmission coefficients, 



we can then use it as a tuning spin-transition parameter. 
In the panel a) of fig. [3] we have the transmission coef- 
ficients J^-j- and T^, while in the panel b) we have the 
transmission coefficients Tj,^ and T-j-^. 

It is easy to see the influence of the rotation opera- 
tor e-^Ebify+v^Cs j not on iy on we n known shifting 

between T-j-^ and 2+^, but also on the spin transition 
— >l processes. The shifting phenomena is given by 
—ialEbCs&z influence, while the spin transition depends 
on —ialEbCs&y The general trend of T-j-^ and Tj,^ is os- 
cillatory. The minimum values of Tj,^ are around zero and 
the minimum values of T^^ are not. The transmission co- 
efficients T-j-j, and T^, exhibit the oscillatory behavior of 
the two uncoupled-spin transmission coefficients. These 
maximum and minimum reflect the passage of flux from 
one spin state to another. The origin of these transitions 
is the precession term of e _la/B i>( cr! ' + ' T ^ s , which stimu- 
lates a mixing between the propagating physical channels. 
The minimum values of T-j-j, and Tj.t are close to zero too. 



Here, we show that by varying £b, the transmission can 
be modulated more efficiently: we can flip the spin of the 
incident electrons, or block the spin-down completely, and 
thus establish a spin transistor. 

An important question is how robust the results are if we 
change the value of the gate electric field £* by £* ± S£ g , 
where = ±0.06. In fig. @] we plot the transmission 
coefficients as function of the bias electric field: in panel a) 
T-tyj-, Tj_4 and are plotted for = +0.06, in panel b) 
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Tff, X 1 ^ and are plotted for -rrf = —0.06. As shown 

g 

there, the transmission curves are not always perfect: the 
minimum values of Tj,^ are not zero and their position 
change. However, the general trends discussed in fig. [3] 
are preserved and, especially, one can modulate them by 
varying the bias electric field. 

All of the results presented so far are valid when only 
a single Fourier s mode propagates in the base. If more 
modes are allowed to mix, when a 7^ j3 and with a mag- 
netic field, the coefficient transmissions pattern become 
more complex but it is still possible to have an analyti- 
cal solution by using the Sylvester Theorem for matrix- 
valued functions for solving the Pauli eq. Details will be 
given elsewhere. In ref. [3] was discussed the multichan- 
nel transmission coefficients for a FET system using this 
approach where the SOI is not relevant. 

In summary, we combined the spin precession in the 
base of a spin-FET, due to the induced spin-orbit coupling, 
with the analytical solution of the Pauli eq., and applied 
it to calculate the transfer matrix at the base. We showed 
that we can select the spin of the outgoing electrons to 
be the same as or opposite to that of the injected spin 
polarized electrons. More important, we can have a nearly 
binary square-wave transmission spin- valve effect for the 
spin-up orientation. 

* * * 
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